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ABSTRACT 

High angular resolution VLBI observations of Centaurus A have been un- 
dertaken that allow access to a wide field-of-view, encompassing both the well- 
studied pc-scale jet and the inner part of the kpc-scale jet. The VLBI observa- 
tions have detected compact regions of synchrotron emission in the kpc-scale jet 
that coincide with three stationary features identified from previous VLA mon- 
itoring observations. Each of these stationary features is associated with strong 
localised X-ray emission. The VLBI results strengthen arguments made by pre- 
vious authors suggesting that the stationary features may be the result of stellar 
objects or gas clouds traversing the jet flow, intercepting the jet and causing 
strong shocks. The VLBI data show that the most strongly shocked regions in 
these features are resolved but have extents no larger than a few pc, reducing 
the required mass of the typical intercepting object by a factor of ~10 relative 
to previous estimates, making explanations based on high mass loss stars or low 
density gas clouds more plausible. 

Subject headings: Galaxies: active. Galaxies: individual (Centaurus A), Galaxies: 
jets, Radio continuum: galaxies 



- 3 - 



Introduction 



Centaurus A, the Fanaroff-Riley type I radio galaxy liosted in the nearby galaxy NGC 
5128, is the closest classical radio galaxy to us and has been the target of a number of 
high angular resolution investigations with a variet y of in s trume nts operating across the 
electromagnetic s pectrum. For a gen eral review, see 



(3.84±0.35 Mpc; 



Rejkuba. M. 



Israeli ( 119981 ) . As a nearby radio galaxy 



(120041 )). high resolution studies of the active nucleus of the 
galaxy are of spe cial interest, in order t o probe as closely as possible to the 5.5 x 10'' solar 



mass black hole (ICappellari et al. 



20091 ). with the aim of gaining insight into the processes 



giving rise to relativistic jets (1 mas ~ 0.02 pc at 3.84 Mpc). 

The active nucleus of Centaurus A has been the subject of comprehensive very 
long baseline interferometry (VLBI) observations, the highest angular resolution direct 



the sub-parsec-scale jet origina 


ting at the galaxy nucleus were explored (Horiuc 


li et al. 


2006: Tineav & Mi 


LirDhv 2001: 


Tineav. Preston & i 


auncev 2001: Fuiisawa et al. 


2000: 


Tineav et a 


1998; 


van Lane;evelde & Beaslev 1998: 


Kellermann. Zensus & Cohen 1997: 


Jones et al. 


1996). 


These probes of the relativistic jet material close to its source have 



yielded information on the speed of flow of the material, the bipolar nature of the flow 
(a jet and counterjet), the size of the "core" of the radio source relative to the expected 
Schwarschild radius of the black hole, and the nature of the ionised and molecular 
environment of the active nucleus. 

In recent years the VLBI technique has also been used to help determine the nature 
of the jets originating from active nuclei on the vastly larger kiloparsec scale. VLBI can 
be used to detect compact radio emission from shocks in radio galaxy iets a f ter th ey have 



propagated large distances from the active nucleus. Recently, 



Tingay et al. 



(l2008h have 



used VLBI observations to detect the compact shocked structures in the termination hot 
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spot of the closest Fanaroff-Riley type II radio galaxy, Pictor A, providing the highest 
resolution view of a radio galaxy hot spot to date and providing new empirical evidence 
that leads to a mor e natural interpretat ion of the X-rays from the Pictor A hot spot than 
previously possible. 



Tingay et al. 



(120081 ) also review other examples of VLBI observations 



of radio galaxy termination hot spots. 

In between the active nucleus and the termination hot spots in radio galaxies, the jet 
propagates large distances and it is supposed that brightness enhancements in the jets seen 
at arcsecond-scale resolution with instruments like the Very Large Array (VLA) are locally 
shocked regions, due to internal instabilities in the flow of material or due to interactions 
with ambient material external to the jet. These enhancements are typically of much lower 
brightness temperature than the radio emission observed from either the active nucleus or 
the termination hot spots. However, a sensitive VLBI search for compact structure in these 
shocked regions may be of great use in understanding the propagation of the jets and the 
processes by which they are shocked. 

Here we describe such a search for compact structures in the kiloparsec-scale jet of 
Centaurus A. The jet on these scales has been y e ry we ll studied with the VLA over a long 



period of time. Most recently, iHardcastle et al. 



(I2OO3I ) combined a multi-epoch study of 



the kiloparsec-scale jet of Centaurus A with the VLA and new Chandra X-ray data. The 
proximity of Centaurus A is crucial in this study, as Centaurus A is the only radio galaxy 
for which the Chandra angular resolution is comparable to the energy-loss travel distance 
of the X-ray emitting electrons. Centaurus A is therefore an importa nt case study for 
understanding jets from active galaxies and the 



Hardcastle et al. 



(120031 ) results reveal much 



about the dynamics of the jet on the kiloparsec-scale, based on radio and X-ray observations 
of the shocked regions of the jet. VLBI observations can add further information, by 
probing the conditions within the shocked regions at high resolution. 
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In section 2 we present our VLBI observational data and imaging results. In section 3 



we discuss the implications of our detection of compact structu res in the shocked 
the kiloparsec-scale jet, relating them to the previous results of iHardcastle et al. 



region s of 



fl2003h 



2. Observations and data analysis 

2.1. LB A Observations 

A VLBI observation of Centaurus A (PKS B1322-428) was made on 21 June, 2007 
using a number of the Australian Long Baseline Array (LBA) telescopes: the 70 m NASA 
Deep Space Network (DSN) antenna at Tidbinbilla (5 hr track, RR polarisation only); 
the 64 m antenna of the Australia Telescope National Facility (ATNF) near Parkes (10 
hr track); 5 x 22 m antennas of the ATNF Australia Telescope Compact Array (ATCA) 
near Narrabri used as a phased array (12 hr track); the ATNF Mopra 22 m antenna near 
Coonabarabran (10 hr track); the University of Tasmania's 26 m antenna near Hobart 
(9 hr track); and the University of Tasmania's 30 m antenna near Ceduna (11 hr track). 
The observation utilised hard disk data recorders (Phillips et al. 2009, in preparation); 
allowing dual circular polarisation data across 4 x 16 MHz bands (2268-2332 MHz) to 
be recorded at Parkes, ATCA, and Mopra, and 2 x 16 MHz bands (2268-2300 MHz) at 
the remaining antennas. Observing scans were centred on the nucleus of Centaurus A 
(a = 13'^'25°>275615200; 6 = -43°1'8'.'80500 [J2000]) which doubled as a phase reference 
source. 



2003) 



The data were correlated with the DiFX Software Correlator (iDeller et al. 
using an integration time of 2 seconds and 64 frequency channels across each 16 MHz 
band (channel widths of 0.25 MHz). The ATCA primary beam limits the field of view of 
this observation to a half width at half maximum (HWHM) of ~ 47". At the HWHM 
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point, bandwidth smearing (a form o f chromatic aberration) and time- averaging smearing 



( ICotton 



1999 



Bridle &: Schwab 



19991 ) fosses are estimated to be approximately 5% and 



respectively. 



2.2. LBA Data Reduction 

The correlated data were imported into the AlP^ package for initial processing. 
The data were fringe-fit (AlPS task FRING) using a one minute solution interval, finding 
independent solutions for each of the 16 MHz bands. 

The nominal system temperatures applied for each antenna during correlation were 
refined by applying the antenna system temperatures measured throughout the observation. 
Further refinement in the calibration was complicated by t he complex jet structure of 



Cent aurus A. To accou nt for this struct ure a new DIFMAP (jShepherd. Pearson and Taylor 



19941 ) task, cordumpQ (ILenc et al. 



20061 ). was used to transfer all phase and amplitude 
corrections made in DIFMAP during the imaging process to an AIPS compatible solution 
table. The Centaurus A data were averaged in frequency and exported to DIFMAP where 
several iterations of modelling and self-calibration of both phases and amplitudes were 
performed. The resulting contour map of the jet and counter-jet are shown in Figure [H A 
one sigma RMS noise of 0.27 mJy beam~^ is measured in the final image; just over twice 



^The Astronomical Image Processing System (AIPS) was developed and is maintained by 
the National Radio Astronomy Observatory, which is operated by Associated Universities, 
Inc., under co-operative agreement with the National Science Foundation 

^The cordump patch is available for DIFMAP at 

http : / / astronomy . swin . edu . au/$\sim$elenc/Dif mapPatches/J 
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the theoretical thermal noise for the observatioiu. The phase and amplitude corrections 
were transferred back to the unaveraged AIPS data set with cordump and the jet model 
components were subtracted from the calibrated data set (AIPS task UVSUB) to simplify 
detection of sources in the wider field. 



To reduce non-coplanar effects (jPerley 



19991 ) ■ the visibilities were {u,v) shifted (AIPS 



task UVFIX) approximately 16" north-east of the nucleus so that the new phase center 



20031 ). Initial 



would coincide with bright knots observed with the VLA (iHardcastle et al. 
imaging of the knot region, using all baselines and uniform weighting, revealed no detections 
above a 6a detection threshold suggesting that sources in this region were either extended, 
weak, or both. A second imaging attempt included only the shortest baselines (between 
Parkes, ATCA and Mopra) and natural weighting to improve sensitivity to larger scale 
structure. The lower resolution image revealed four detections above the 6a threshold 



[Figure [2]) and achieved a one sigma rms noise of 0.25 m. 



produced for display using the KARMA software package (iGooch 



y beam ^. AH images were 



19961). 



Results 



3.1. Identification of Sources 



A list of sources detected above the 6a detection threshold of the short-baseline VLBI 
data is given in Table [H The 6a threshold provides a conservative false detection rate of 
less than one in 3000. Flux density errors of ±10% are expect ed due to uncer tainties in 



the absolute flux density scale for Southern Hemisphere VLBI (IReynolds 



1994h . The total 



flux density and sizes of each source were determined by modelling each component with 



■^Estimated with the ATNF VLBI sensitivity calculator: 
'http : //www. atnf . csiro . au/vlbi/calculator/ 
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elliptical Gaussian components, the peak flux density, integrated flux density, fltted size 
and position angle (PA) of each source is listed in Table [TJ 



Component 


Size 


P.A.^ 


S'2.3 Peak 


5*2.3 Integrated 


Tb 




(mas) 


(deg) 


(mJy beam^^) 


(mJy) 


(K) 


AlAa 


83 X 39 


-73 


3.8 


6.1 


4.3 X 10^ 


AlAb 


109 X 45 


47 


3.1 


9.5 


4.5 X 10^ 


AlC 


260 X 70 


-59 


2.4 


9.8 


1.2 X 10^ 


A2Aa 


106 X 106 


-40 


2.2 


8.0 


1.6 X 10^ 


A2Ab 


320 X 64 


-83 


1.5 


5.6 


6.3 X 10^ 



Table 1: VLBI component parameters. 



"Position angle of major axis measured east of north. 



Figure 2 sho ws the detection o : 



An overlay of the 



Hardcastle et al 



the weak but compact radio sources listed in Table 1. 



(I2OO3I ) image and our VLBI image (Figure 3) shows that 
the weak, compact radio sources coincide with features AlA (resolved into subcomponents 
we designate AlAa and AlAb), AlC and A2A (resolved into subcomponents we designate 
A2Aa and A2Ab). These compact sources are each better than seven times the RMS noise 
level in the VLBI image of 0.25 mJy/beam and are t herefore detected with go od confldence. 



Hardcastle et al.i (i2003n . one other 



On the basis of the VLA flux densities presented in 
component may have been expected to be detected, AlB. AlB was not detected in our 
VLBI data. 
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4. Discussion 



Hardcastle et al.l (120031 ) examine a number of different possibilities for tlie nature of the 
discrete features in the jet of Centaurus A, detected at both radio and X-ray wavelengths. 
We examine the new VLBI imaging results in the light of these possibilities. 



Hardcastle et al.l (120031 ) briefly consider the possibility that the features are jet- 
associated radio supernovae or supernova remnants. Certainly at the distance of Centaurus 
A, the compact objects we have detected with VLBI, at a few mJy each, are typical (in 
terms of flux density) of radio supe rnova remn ants detected in nearby starburst g alaxies 



such as M82 (iMcDonald et al. 



4945 (ILenc fc Tingay 



200l|), NGC 253 (ILenc fc Tingav 



2006 



Tingay 



20041 ) or NGC 



20091 ). While very few supernova remnants seen in other galaxies 



can be resolved into a shell morphology, some clearly are. The fact that the compact 
sources detected in Figure 2 do not exhibit shell morphologies is not necessarily evidence 
against them as supern ova remnants. The evid ence against these sources being supernova 



remnants, as argued by 



Hardcastle et al. 



(120031 ) . is that the X-ray emission associated with 



these radio components is apparently non-thermal. 



Hardcastle et al 



(120031 ) also argue against the possibility that the jet features are 
purely compressions in the synchrotron-emitting fluid, suggesting that the features instead 
represent privileged sites of in situ particle acceleration. These sites may be standing shocks 
in the fluid flow, caused by the jet's interaction with a stationary obje ct of some descrip tion 



high mass loss stars or gas clouds are cited as possible examples. 



Hardcastle et al. 



( 120031 ) suggest that the Al complex of components may represent standing shocks related 



to the transition from the collimated inner jet to the flared extended jet, reserving the 
interpretation of jet interaction with stationary objects for the stationary components 
(stationary simply meaning that no proper motion has been detected) downstream of 
the Al complex (for example component A2A). However, it is apparent from Figure 1 in 
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Hardcastle et al.l (120031 ) that the inner jet commences flaring well before the jet reaches 
the Al complex. Why components AlA and AlC should therefore be stationary shocks at 
that particular point in the jet is not clear. Under this interpretation it is also not obvious 
why component AlB should have a large proper motion as part of the complex of shocks 
at the jet transition point. A common interpretation for all of the stationary components 
in the jet therefore remains a possibility, that all stationary radio features are related to 
interactions with objects that have intercepted the jet flow. 

From our VLBl observations we find tha t we detect compact rad io emission from the 



Hardcastle et al. 



(I2OO3I ) Chandra data, these 



VIA features AlA, AlC and A2A. From the 
features correspond to AXIA, AXlC and AX2. These are all significant X-ray sources 
and all have radio counterparts in the VLA image that are stationary. The VLBI data 
now show that these three features have radio emission compact on parsec scales. Simply 
on the basis of the VLA radio flux densities as a predictor of compact emission on VLBl 
scales, one other compact radio component may have been expected to be detected with 
VLBI, along with AlA, AlC and A2A, that is AlB. AlB was not detected with VLBI, 
although it is of comparable strength and size to the other three on VLA scales. AlB is 
the only featu re of the four tha t has a large apparent motion away from the nucleus. The 



proposition of 



Hardcastle et al 



( I2OO3I ) that the stationary features in the Centaurus A jet 
represent stationary shocks is strengthened by our VLBI results that show compact radio 
structures associated with these stationary features. 

The sizes of the compact radio sources detected with VLBI are all less than ^6 pc 



Hardcastle et al. 



(typic al size of ~2 - 3 pc), significantly less than the 10 pc assumed by 
(120031 ) when calculating the properties of the objects that could possibly explain the 
presence of the stationary shocks. The sizes of the sub-structures seen in the AlAa and 



AlAb components at VLBI resolution are consistent with the sizes inferred from ~10% 
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flux density variability over a decade detected in the VLA data (iHardcastle et al.ll2003l ). 
We note that due to sensitivity and surface brightness sensitivity limitations with the 
VLBl data, weaker and/or more extended emission could be associated with the compact 
components, the structure in the VLBI images perhaps representing the core of the shocked 
region. However, the primary interactions in the shocked regions are likely to be of order 
only a few pc in extent. Taking the VLBI sizes for the shocked regions as indicative of the 
sizes of the interaction regions in the jet implies that the constraint on the required mass in 
the obstacles is reduced by a factor of ~10, compared to estim ates based on the comp onent 



sizes as measured with the VLA (see equation in section 5.3 of 



Hardcastle et al. 



(120031) that 



has the mass of the obstacle proportional to the square of the radius of the obstacle). 

This reduction in the required mass of the obstacle, from at least a few solar masses to 
su bstantially less than one solar mass, relaxes some of the extreme characteristics required 



by 



Hardcastle et al 



(120031 ) when discussing obstacles based on high mass loss stars and low 



density gas clouds 

The possibility that the stationary shocked regions could be due to gas clouds 
intercepting the jet is a scenari o that has been st udied by a number of authors using 



numerical simulations. Recently 



Jeyakumarl (120091 ) reported hydrodynamical simulations of 



jets interacting with clouds of various sizes that move through the jet at various speeds, 
showing the amount of disruption caused. Judged qualitatively against the images showing 
the kpc-scale morphology of the Centaurus A jets, one could conclude that the stationary 
features in the kpc-scale jet are due to relatively fast moving clouds (~10,000 km s~^) with 
diameters only a small fraction of the jet width. The VLBI results can be taken to support 
such a conclusion - a cloud diameter of a few pc covering approximately 5% of the jet width. 

A comprehensi v e rec ent study of jet interactions with stationary clouds by 
Choi. Witta &: Ryul (120071 ) shows that compact synchrotron emission can be generated 
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from these types of interactions. Again, qualitatively, the compact synchrotron-emitting 
structures seen in the numerical simulations appear similar to those observed with VLBI, 
in particular components AlAa/b and A2Aa/b in Figure 2 that show a bright compact 
feature and a weaker and more extended trailing feature. 

5. Conclusions 

We have made the first VLBI images of pc-scale structures in the kpc-scale jet of 
Centuarus A, detecting compact regions of synchrotron emission. All detected pc-scale 
components coincide with the locations of jet components seen to be stationary from 
previous VLA monitoring of other authors. Other features in the kpc-scale jet have 
been shown to have large proper motions, from the VLA monitoring. None of these 
moving components have been detected with VLBI. These results support suggestions 
that the stationary components in the jet may be due to moving stellar objects or gas 
clouds intersecting the jet flow and causing strong shocks. Size estimates for the shocked 
regions suggest extents of order a few pc, smaller than previously assumed and consistent 
with observed flux density variability, leading to approximately an order of magnitude 
reduction in the required mass for the obstacles. This reduction relaxes some of the extreme 
constraints on the inferred properties of high mass loss stars and low density gas clouds 
suggested as providing plausible obstacles. A repeat of these VLBI observations at a later 
date will put limits on proper motions for the detected pc-scale components in the kpc-scale 
jet of Centaurus A, allowing further study of their properties. 

This research has made use of NASA's Astrophysics Data System. This research has 
made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under contract with the National 
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Fig. 1. — Uniformly-weighted 2.3 GHz LBA image of the jet and counterjet in Cen A. 
The peak surface brightness is 1.46 Jy beam~^ and the la rms image noise is 0.27 mJy 
beam~^. Contours are drawn at ±1, ±2, ±4, ±8, • • • times the 6a rms noise. The beam size 
is 10.8 X 10.5 mas at a position angle of 37°. 
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Fig. 2. — Naturally-weighted 2.3 GHz LBA image of knot features in Cen A. The peak 
surface brightness is 3.8 mJy beam~^ and the la rms image noise is 0.25 mJy beam~^. 
Contours are drawn at ±2^^, ±2^, ±2^, ±2i, • • • times the 3a rms noise. The beam size is 
150 X 37 mas at a position angle of —62°. 
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Fig. 3. — Uniformly-weighted 8.4 GHz VLA contour map of Centaurus A (grey contours) 
overlaid with a naturally- weighted 2.3 GHz LB A contour map (black contours) to show the 
locations of the compact sources detected within the jet and knot features. 



